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Abstract. Adopting a quantum information perspective, we analyse the correlations in the 

thermal light beams used to demonstrate the Hanbury Brown and Twiss effect. We find that the 
00 total correlations measured by the Renyi mutual information match the normalised intensity 

correlations in the regime of low source intensity. Genuine quantum correlations in the form 
,——1 of discord are relevant in such regime but get washed out with increasing source intensity. This 

provides a new angle on the issue about the nature — quantum versus classical — of the effect. 

^ 1. Introduction 
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Figure 1 : Diagrammatic representation of a Hanbury Brown and Twiss set-up. Light from a large source 
propagates into the far field whereupon it is split on a balanced beamsplitter. We denote by a the mode operators 
^ of the source, by Uvac the mode operators of the vacuum field entering the other arm of the beamsplitter, and by 

bi2 the outcoming mode operators in the two arms after the beamsphtter. When one of the detectors is scanned 
in the transverse plane, the normalised intensity correlations of the modes bi 2 exhibit an oscillatory behaviour 
as a function of the transverse coordinate x of the scanning detector. We find that these intensity coiTelations 
are matched by the total quadrature correlations (as measured by Renyi-2 mutual information) in the bipartite 
system formed by the modes bi 2 in the limit of low source photon-count. In that regime, a relevant portion of 
such total correlations is constituted by genuinely quantum ones as measured by discord. 



The Hanbury Brown and Twiss (HBT) effect was demonstrated in 1956 [IJ and had 
far-reaching impact, even inspiring some of Glauber's seminal work on quantum optics f2^. 
It was the first demonstration of the second-order correlation of visible light and caused 
significant debate; some held the belief that it ran counter to the photon theory of light and the 
experimental results were therefore incorrect [I21lll[5]|. After a series of papers by Hanbury 
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Brown and Twiss, the controversy was largely abated [[6ll21. However, even in the past few 
years, papers have been published which argue against the established interpretation (Hlllll. 

This revived interest comes about from an analogy between the HBT effect and 
the phenomenon of ghost imaging: simply put, both phenomena depend upon intensity 
correlations. Due to questions surrounding the classical interpretations of ghost imaging, 
similar suspicions were raised upon the HBT effect by Scarcelli et al. who suggested that its 
explanation relies on non-local two photon interference [8J, although this claim has proved 
controversial IfTOl . 

The HBT effect can be demonstrated by a set-up as shown in figure 1 . Light from a distant 
source impinges upon a balanced beamsplitter and then it propagates down two separate, yet 
identical, arms of the configuration. Subsequently, each beam falls upon a point-like detector. 
When one detector is kept fixed and the other detector is scanned in the transverse plane, we 
observe a correlation of the intensity in what is known as the HBT effect. 

Prior to the 1956 HBT paper, such effects had been demonstrated in radio astronomy and 
were readily accepted as permitting a simple, classical explanation [jT2l. This is because the 
experiments had all been performed in the regime of large photon-number per mode, in which 
the classical theory of light easily holds sway [031 and the shot noise does not contribute a 
large extent. In this case, we can regard the beamsplitter as creating two identical duplicates 
of the intensity profile of light. Then by the van Cittert-Zernike theorem it is possible to 
classically calculate the correlation function as we scan one detector across the beam [[TT]|. 
Essentially, the light propagates in such a manner that the classical intensity fluctuations are 
appreciably correlated at the detection plane over a finite area. 

However, the 1956 experiment involved measuring the correlations of optical photons 
from a mercury arc lamp, in which each mode contains few photons. In this case, it becomes 
more difficult to understand how it can be that some sort of 'classical' behaviour emerges even 
when few photons impinge upon the detectors. However, it was proven by Purcell [4], that 
the experimental results of HBT could be explained by using a classical light field as long as 
the electrons in the detector were given a quantum treatment. 

The success of this 'semi-classical' model later came to be understood through the 
generalised formalism of quasi-probability distributions [illil . In 1963, during the birth of 
quantum optics in its current form, Sudarshan proved that any experiments involving light 
which is characterised by a positive, well-behaved ^-representation can also be fully explained 
by semi-classical methods [fT4ll . Thermal light - such as that present in the HBT experiment 
- satisfies this constraint. 

The notion of classicality which emerges from the ^-representation (P-classicality), 
however, is by no means the only notion available. In quantum information, a popular and 
alternative definition of classicality for bipartite states arises from the concept of 'discord', 
a quantity which is defined as the difference between two classically identical expressions 
for mutual information [iTTl [T8l . The discord can be considered to represent the degree of 
quantum correlations between each part of the bipartition and when the discord is equal to 
zero, this defines an indication of classicality (D-classicahty) of the correlations in the state. 
It has recently been shown that both notions of classicality are maximally inequivalent [[T5l . 
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In this paper we follow the example of a recent treatment of ghost imaging in terms of 
discord [fT6l . Since the HBT effect is entirely contingent upon the correlations between each 
arm, we believe that the very general informational approach to correlations as defined by 
the mutual information (total correlations) and discord (quantum correlations) provides an 
apt foundation for studying (non)classicality within the scheme. In particular, we exploit 
a powerful new variant of mutual information (and its associated classical and quantum 
correlations) defined using the Renyi-2 entropy for Gaussian states [19|. This definition of 
mutual information is intimately related to phase space sampling by homodyne detections 
which serves as a strong interpretational aid as well as harbouring numerous other desirable 
properties [fTPTl . 

We establish that, despite the availability of a semi-classical light model, the quantum 
correlations are indeed significant in the regime of low photon-count per mode in which the 
original controversy was conceived. Thus, whilst P-classicality holds, we find an undeniable 
quantum signature as revealed by a relevant portion of quantum correlations (discord). This is 
somehow reflected in the physical model proposed by Scarcelli et al. [8] in which non-local 
correlations feature in the imaging. We further find that in the limit of a very low photon-count 
per mode, the mutual information tends to coincide analytically with the normalised intensity 
correlations. 

2. Preliminaries 

2.1. Renyi entropy information measures for Gaussian states 

A detailed account of the correlations and information measures in terms of Renyi-2 entropy 
can be found in [|19il . To summarise, the definitions reflect the usual formulation of mutual 
information, classical correlations and quantum discord [flTlfTSlI but with the replacement of 
the von Neumann entropy with the Renyi-2 entropy. 
A Renyi-a entropy is defined by 

S, = {l-arHntr(p") (1) 

The von Neumann entropy which is used as a conventional measure of quantum information 
can be obtained by taking the limit a ^ I, whereas for the purposes of this paper we use the 
Renyi-2 entropy, which is defined by a = 2, and simply reduces to minus the logarithm of the 
purity of a quantum state. 

For this choice, and given a bipartite Gaussian state Pab with marginal states Pa (Pb) for 
subsystems A (B), the Renyi-2 mutual information, classical correlations and quantum discord 
are defined respectively by 

I(A:B) = S2(Pa) + S2(Pb)-S 2(Pab) (2) 

J(A\B) =S2(Pa)-M'H2(A\Bu,) (3) 

n/ 

D(A\B) =S2(PB)-S2(PAB) + inf^2iA\Bn,) (4) 

Here, we have introduced a quantum conditional entropy 'K2(^|5n,) = Zr PiS2{pA\d vvhich 
quantifies the entropy on A after a measurement {11,} has been performed on B. Herein 
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lies the subtlety which distinguishes the quantum from classical case. In the classical case, 
the conditional entropy factors out, and the discord evaluates to 0. In the quantum case, 
however, we need to consider the measurement performed when calculating the conditional 
entropy. Since a measurement can disturb the system, the means we cannot fully extract 
all the information of a quantum system. Taking the least disturbing measurement, and 
restricting this minimisation to Gaussian measurements, we can calculate the Gaussian 
classical correlations and quantum discord V as defined in the above equations ([3jj4]) 

mm. 

We remind that the mutual information I defined in terms of the Renyi-2 entropy as 
in equation (|2]), of a particular Gaussian state p^e, corresponds to the Shannon continuous 
mutual information of its Wigner distribution. More generally, it stems from the satisfaction 
of the strong subadditivity inequality by the Renyi-2 entropy (on Gaussian states), that it can 
be used to reformulate the core of quantum information theory and to define valid correlation 
measures in the Gaussian setting [19]. 



2.2. Covariance matrix description of the HBT effect 

In order to compute the discord and other correlations for Gaussian states, it is necessary to 
produce the covariance matrix of the relevant quadrature operators fSOl. Such a matrix fully 
characterises the state up to local displacements. For a general two-mode Gaussian state Pab, 
any covariance matrix (Tab can be written (by means of local unitaries) in the standard form, 

^ a c ^ 



o-AB = 



a d 
c b 
d b 



(5) 



For the HBT setting under consideration (see Fig.[T]), the entries are easily calculated with 
some simplifying assumptions: firstly, we choose a large disk-like monochromatic source of 
area A, emitting spatially incoherent light. From this, we take each mode a(q) (where q is 
the transverse component of the wave- vector) to be uncorrelated with every other mode in the 
source plane, and we consider all modes to be zero-mean modes with identical photon-number 
expectation value n, such that {a\q)a{q')) = n6{q - q'), these are common assumptions 
for intensity interferometric experiments with thermal light [|8l|22]|. This implies that in the 
source plane {a'^(x2)a(xi)) = J e''^'^^e'^''^^{a'^(q')a(q))dqdq' = h6{xi - X2), where Xi represent 
transverse position vectors. Furthermore, we take the beamsplitter to be balanced (50:50) 
such that the outcoming mode operators for each branch are 

1 

bi,2(^) = -^{a{x) + Kacix))- (6) 

Now we can calculate the covariance between all four mode operators of relevance bt and b\, 
with / = 1,2 (after which we can rotate into the quadrature basis). To do this, we also need 
to take into account the propagation of the incoming modes a{x) through space. We consider 
the Fraunhofer propagator g{x\,x) = g 'i*^' which relates the modes in the source plane to 
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the observation plane by a(x) oc f^^^^^^^ g{xi,x)a(xi)dxi. As an example, we can then calculate 
and report here the correlation function {a\x)a(x')). 



(a^(x)a(x)) oc 



g{xi,x)g*{x2, x')(a(xi)a(x2))d^xid^X2 



n 

J A J A 

//' 

J A J A 



,iz(x2-x' -xi-x)- 



n6{x\ - X2)d x\d X2 



2j{hA\x-x'\) (kA 
= n — = n Jinc 

^-^\x-x'\ 



X - X 



(7) 



where J refers to a Bessel function of the first kind and we have defined Jinc(y) = 
Continuing along these lines and setting x' = we obtain a covariance matrix auix) in 
standard form as in equation ([5]), where 



a = b = I + 2n, and 
kA 

c = d = 2n Jinc | — \x\ 



(8) 



2.3. Measures of correlations 

From the previous analysis, and using the results of [fT9l , we can obtain the explicit 
expressions for the Renyi-2 mutual information, classical correlations and discord as functions 
of X, 

(9) 



I{A : 5) = In 
J{A\B) =ln 
!D{A\B) =ln 



a^ - c^ 



a + a 
a^ + a - c^ 



9 T 9 

a + a — ac 



(10) 



(11) 



a^ + a^ - ac^ - c^ 

We find it relevant at this point to recall also the formula for the normalised intensity 
correlation function in our scheme. Defining the intensity in each beam / = 1,2 as 
Ii(x) = b]{x)bi{x), we have the well established result IfTTll 

2 



{h{x')h{x)) 

{h{x')){l2{x)) 



- 1 



,kA 
Jmc I — \x - X 

z 



(12) 



3. Results 



To begin with, it is interesting to see how the classical and quantum correlations behave as 
we scan the detector for fixed values of h. In all the following results, we have normalised 
the total correlations to a maximum value of 1 by plotting ; and similarly for quantum and 
classical ones. For large h as in figure [2[^a), it is apparent that the quantum correlations decay 
faster than the classical ones and also contribute a much smaller part of the total correlations 



Unveiling the Hanbury Brown and Twiss effect through Renyi entropy correlations 



6 



even at their peak. However, for very small values of h as in figure |2jb) this is no longer the 
case: we see that both classical and quantum correlations degrade at the same rate and are 
approximately equal. This result stems from the fact that the purity of a Gaussian state relates 
only to the number of thermal photons [,23 J present and for small photons counts, we obtain 
high purity. Recall that for pure states, D = = ^I. 

To better understand the behaviour of the correlations as we vary h, it is useful to plot 
changes against n at fixed x, see Fig. ^c). It is apparent that the total correlations are 
composed of a much greater portion of classical correlations than quantum correlations at 
high photon-count per mode. This gels closely with our intuition: in the regime of high 
brightness, the quantum component of correlations becomes very small. It is thus easy to 
see why in certain regimes, the effect is simple to explain classically: it appears that high 
photon counts wash-out quantumness, increasing the mixedness and simultaneously quashing 
the presence of discord. 

Returning to the regime of high purity (and low photon-count), we find a notable 
result. The total correlations as quantified by I match almost exactly the normalised 
intensity correlations of equation ( [T2] ). Notice once again that classical and quantum 
correlations (discord) both contribute in equal halves to such total correlations, and thus to the 
manifestation of the HBT effect itself in this regime. Very remarkably, for any form of cross- 
correlation of narrowband, thermal light, the correspondence between mutual information 
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Figure 2: Plots of the normalised total (I, solid black Hne), classical (JT, dashed red line), and quantum (£), 
dotted blue line) correlations as a function of x for (a) h - 10 and (b) h - 0.01; the thin green hne represents 
the normalised intensity correlations, which practically coincides with the mutual information in (b). Panel (c) 
depicts the correlations as a function of n with fixed x - Q; the inset shows the same correlations normalised by 
the mutual information, to better highlight the quantum and classical contributions to the total correlations. In 
all the plots kA/z = 1000. 
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and normalised intensity correlations holds analytically in the weak-light regime (up to the 
third order in h), as proven in the Appendix. This provides an intriguing connection and a 
convenient short-cut for calculating and/or measuring experimentally these correlations. 

4. Discussion and conclusions 

We have provided an example of how the Renyi-2 entropy [[T9B can serve as a powerful tool 
for analysing the behaviour of correlations in the Gaussian world. The presented results, 
visually summarised in Fig. [2} encapsulate not only information about the correlations but also 
provide, at a mere glance, information about the purity. They also reveal, significantly, that in 
low-illumination conditions the normalised mutual information can be properly approximated 
by the normalised intensity correlation. A somehow similar result has been recently found in 
the context of thermal light ghost imaging, whereby a properly defined coarse-grained mutual 
information has been found to match the signal-to-noise ratio of the scheme in the limit of 
low illumination [16l. 

Applying the tools of quantum information to the HBT eff"ect clarifies the source of 
the disagreement in its interpretation and unveils a quantifier of quantumness which behaves 
particularly intuitively. We have shown that in the limit of a high photon-count per mode the 
correlations are mostly classical, which matches observations made as far back as the 1950s 
that the classical 'wave-like' features of light emerge in this regime 02411 . 

We find that for very low photon-counts, however, the quantum correlations are no 
longer negligible and it was in this regime where the effect originally caused ripples of 
controversy. This provides a significant example of an insufficiency of following a single 
definition of classicality. Whilst Scarcelli et al. interpreted their finding of non-local two- 
photon interference as contradicting previous semi-classical interpretations [8J, others have 
shown this to be incorrect [|25l . The actuality of the situation, and the source of confusion, is 
that the P-classicality criterion is simply unable to resolve the quantumness in the system. 

On the other hand the present discord analysis does support the finding that there is 
a significant presence of quantum effects in the light correlations at low illuminations. As 
previously mentioned, the strong discrepancy between both definitions of classicality has 
been recently delineated on rigorous footings [IT5l|. It is sufficient for our purposes though 
to note that the P-classicality is equivalent to separability for bipartite (squeezed) thermal 
states ||26l while discord captures all quantum correlations beyond entanglement (in fact, any 
two correlated Gaussian states will have non-zero discord [|20ll ). 

In summary, we have utilised a definition of 'quantumness' of correlations based on 
Renyi-2 Gaussian discord [[T9l l20l . which is found to vary continuously from its largest 
relative value in the low-illumination, high-purity regime, shrinking as the illumination is 
increased. We believe that this is a naturally relevant way of looking at classicality for such 
interferometry experiments since the P-classicality reveals very little information and holds 
even in the seemingly non-classical regime of very few photons. 
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Appendix A. 

Here we prove that for two-mode Gaussian state with covariance matrix of the form (|5]) with 

a = b = I + 2h and c = d = 2n fix), we have when n «; 1 that « f(x) = i | = 

^ j(0) J ^ ^ \[(iix)){i{x'))y'^ J 

(h(x'))(h(x)) ~ ^' rightmost equality is well known to hold for thermal light [fTTI . It is easy 

(Q]' (x)b2(x' )) \^ ^ - 

[</(x))(/(v')>]'^^ I = f{x) arises naturally from the fact that {b\{x)b2{x'')) = nf{x) 

and {b.{x)bi{x)) = n where i = 1,2. In order to show that « f(x), it is useful first 

to recall that in standard form we can write I{x) = In -2z~na- write g(x) 



then we can say = |^|^ = log^^Q^ ^(.x:). This then reduces the problem of showing that 
fgj « fix) to showing that g(x) « ^(0)^^''^ The Taylor expansions of (^i+2S''-(£fix))^ ^'^'^ °^ 
[ (i+2«")H2n)^ -^^'^'^^ match up to third order in h (where we have imposed the condition /(O) = 1 
as a natural consequence of the beamsplitter transformation). This provides a proof that in 
low-light source conditions (i.e. small h), the normalised intensity correlations are almost 
exactly equal to the normalised Renyi-2 mutual information in our setting. 



